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Abstract 

We obtain the free energies and critical exponents of models associated with elliptic 
solutions of the star-triangle relation and reflection equation. The models considered are 
related to the affine Lie algebras B 1 ^ , Cn and Dn^ ■ The bulk and surface specific 

heat exponents are seen to satisfy the scaling relation 2a s — a& + 2. It follows from scaling 
relations that in regime III the correlation length exponent v is given by v = (I + g)/2g, 
where I is the level and g is the dual Coxeter number. In regime II we find v = (I + g)/2l. 

KEYWORDS: Exactly solved models, inversion relations, free energy, critical exponents, scaling 
relations. 



The formulation of the boundary version of the Yang-Baxter equation has provided a sys- 
tematic framework for the investigation of integrable models with a boundary [1-3]. Recent at- 
tention has turned to exploiting boundary integrability to derive off-critical surface phenomena, 
such as the surface magnetization p[, the surface free energy and related critical exponents 
[6-11]. For face models, integrability in the bulk is assured by solutions of the star-triangle 
relation (STR) g§ 
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Integrability at a boundary requires the additional relation, 
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which is the face formulation of the reflection equation (RE) |i~3| , |14| , |6|, [T^, 16]. 

Solutions to the RE, defining boundary weights, have recently been found for the An , An 
and X n = Bn^ , Cn^ , models [ 17] , for which the solutions of the STR have been known for 
some time [18, 19|. Our interest here lies in the critical behaviour of the bulk and surface free 
energies of the elliptic face models associated with the algebras and X& . These models 
include some well known models as special cases. For example, the Andrews-Baxter-Forrester 



(ABF) model f2TJ is related to A\ 
Two inversion relations, 
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are satisfied by the bulk weights of the models under consideration. On the other hand, the 
diagonal solutions of the RE found in [|l7j fulfill the boundary crossing relation 




2u + A \K 
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u + A = q s {u)K b 



u 
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The G a are crossing factors jl£|. The crossing parameter is given by A = —tg/2 where the 
parameters t,g are given in Table 1. The key ingredients are the inversion function g(u) and 
the boundary crossing function g s (u). 

The inverson relation method |l^] has recently been applied to a number of models to obtain 
the off-critical surface free energy [7-11]. The unitarity relation 
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for the transfer matrix eigenvalues T(u) follow from the crossing unitarity relation and dis- 



regarding finite-size corrections. Define Tb(u) 



and T s (u) 



then the bulk and the 



log K s 



surface free energies per site can be defined by fb(u) = — logKb(u) and f s (u) 
respectively. 

The restricted-solid-on-solid (RSOS) models follow in a natural way from the unrestricted 
models that we have discussed so far. One introduces a positive integer / and sets L as specified 
in Table 1. Local state a in the Boltzmann weights is taken as a level / dominant integral weight 
of A^ and X$t\ In Table 1 we have also listed the levels under consideration. We do not treat 
/ = 1 for the Bn\Dn^ RSOS models as they are then completely frozen. 



We begin with the ^ model, for which 
[1 + «][!-«] 



and 
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Here we define 
[u] = [u,p] 
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K b (u)K b (-l + u) = " J ^ " J (11) 



k s (u)k s (-1 + u) = ^ ( 12 ) 



where 

oo 

■&1 (u, p) = 2p 1/s sin u JJ (1 - 2p n cos 2u + p 2n ) (1 - p n ) (10) 

n=l 

is a standard elliptic theta- function of nome p = e 2mr . We first consider the model in regime 
III (-1 < u < and < p < 1). 

After taking a convenient normalisation in (^) we have 

[1 + - u] 
[IF 

for the bulk and 

[2 + 2u] [2 - 2u) 
W 

for the surface. To proceed, we introduce the new variables 

x = e -4- 2 A/ £) w = e -A^u/e^ q = e -to*L/e } (13) 

where A = — 1 for the model. The conjugate modulus transformation of the theta-function, 
M™/L,p) ~E{w,q), (14) 
is also required, where p = e~ e l L and 

oo 

E(z, y) = J] (l - y n ^)(i - y n z~ 1 )(i - y n ). (15) 

n=l 

We suppose that n b (w) is analytic and nonzero in the annulus 1 < w < x and Laurent expand 
\ogK b {w) in powers of w. Then matching coefficients in (O) we obtain 



sinh(^) sinh ( 2 " 2 f+^ ) cosh (^M^l 
nsinh(^^)cosh( 
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for the bulk free energy. In a similar manner, we obtain the surface free energy 



sinh(^) sinh ( 4A f+"> ) cosh 



f s (u)= Y Vtt A 2 f -■ (17) 

„iHc nsinh(^)cosh(^) 

Now consider the A J model in regime II (0 < u < — 1 + L/2 and < p < 1). In this case 
we need to modify for the appropriate analyticity strip, with 

( \ i \ [l + u][l-u] 

K b (u)K b (-U) = ( 18 ) 

Kb (u) Kb (L - 2 - u) = [2 + f u] (19) 

[1J 

for the bulk and 

k s {u)k s {-1 + L/2 + u) = 1 - iL (20) 
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for the surface. We assume that n b {u) and n s (u) are analytic and nonzero in this regime, and 
in a similar manner obtain 

oo sin h(^) Sinh ( ^-3) ) ginh ^ 2^n(L-l- U ) ^ 
fb(u) = - ■ ; 7^T7~T (2^n(L~2) _ 



nsinh( i7r nL ) sinh 
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oo S inh(g^) sinh(^) sinh ( 2n n{l+u) ) 
Jlo nsinh(^)sinh(^^)) ' 1 j 

oo sinh(^) Sinh ^ML-2-2u) ^ cpsh ^ 2^n(L-4) ^ 

n=-oo nsinh(^) cosh I 1 J 

for the bulk and surface free energy. 

Now consider the and models in regime III (A < u < and < p < 1) with A as 
given in Table 1. For these models the inversion and bounday crossing functions are given by 

[X + u][X -u][l + u][l -u] 



Q{u) = 1 " JL 7 X ; 2 J 112 JL J ^ (23) 



^)^ [2A + % A - 2 - ] (24) 



After appropriate normalization, we have 

[-A + u][-X -u][l + u][l - u 

for the bulk and 
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for the surface. Under the appropriate analyticity assumptions we obtain the bulk and surface 
free energies 

oo sin h(^) Sinh ( 2 - 2 "^") ) cosh ( *r 2 n(L+A-l) ) cosh 

A(«) = -2j^ ! ns - nh( 2^ )cosh( 2|^ ! > ( 27 ) 

oo sinh (^) sinh ( ^M^u) \ cQgh / 2 ^n(L+3A-l) \ h (2^+1)^ 



nsinh(^^)cosh(^M) 



Now consider the Bn \ Cn^ and models in regime II (0 < u < X + L/2 and < p < 1) 
Similar to the model the inversion relations are modified to 



\-X + u\\-X-u]\l + u\\l-u\ 

K b {u)K b {-u) = [ - JL r _ A]2[ J { ]2 iA i , (29) 

. , n] ) \ [u][-2X + u][l-X + u][-l-X + v : 

K b (u)K b {L + 2X- U) = r ^2^2 ' ( 30 ) 



K („\ K (\ + T/9 + ,A [~2A + 2u][-2X - 2u][l - X + 2u}[l - X - 2u] 

k s (u)k s (X + L/2 + u) = [-2A] 2 [1 - A] 2 ' ^ ' 
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From these relations we obtain 



oo sinh( 2A^) cpsh ( 2^n(A+l) \ / 2^n(L+A-u) \ ^ / 



27r 2 n(L+2A-l) 

e 



nsinh( 2^nL )sinh ( 2^n(L+2A) \ 
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n=-oo nsinh C 1 ^ ) cosh f ^ ' 



for the bulk and surface free energies. Note that in this case we also obtain the free energies of 
the level / ^n-i model, which corresponds to the level n model on changing the signs of 
u and A, as follows from level-rank duality [19|. 

We are particularly interested in the critical behavior of these models as the nome p = 
e -e/L _^ q j n eacri case ^he singular term in the free energies is obtained by making use of the 
Poisson summation formula 



For the model in regime III we have 
| p 2 ~ ab logp for L = 2m 

Jb ~ \ r j . , ( 34 ) 

I nsc tor L = 2m + 1 

for the bulk, where nsc denotes "no singular contribution" and m is some integer. The bulk 
specific heat exponent is given by 

a b = 2 - |. (35) 
For the surface free energy we find 

p 2 ~ as for L = 2m + 1 

fs ~ { p 2 ~ as logp for L = Am (36) 
nsc for L = Am + 2 

where the excess specific heat exponent is given by 

«- = 2 -|- ( 37 ) 
On the other hand, for the model in regime II we have 



, p 2 - ab for L + A 

p 2-a ftlogp for L = 4 



for the bulk, with 

For the surface free energy, 
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with 
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The bulk results (34), (|3^), (23), (pM) have been obtained for the ABF model |2Q] , as have the 
surface results © and © g ||. 



For the Bn^ and Dn models in regime III we have 




for L / -2mA, L 7^ -2mA + 1 
for L = —2mA 
for L = -2mA + 1 
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While for the surface energy, 
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For the l?!^, Cn^ and Dn^ models in regime II we have 



fb ~ < 



nsc 



for LT2X ^ m i and / 2m 2 - 1 



for ZTII = m i and / 2m 2 - 1 

for ^ = - 1 



with exponent 
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For the surface energy 
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2{LT2Ay ^ and ^±g=i ^ 2m 2 - 1 and ^ ^ 2m 3 



2CLT2A) = m i and i rix 1 + 2m 2 - 1 and ^ ^ 2m 3 - 1 (48) 
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with exponent 
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In the above mi, m 2 and m 3 are arbitrary integers. 
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(2) 

In this case the exponents of the level I A^-i model follow under level-rank duality with 
the level n C, model on changing the sign of A. 

The bulk and surface specific heat exponents are seen to satisfy the relation 2a s = 2 + a&. 
More generally, this relation can be inferred directly from a comparison of the singular behaviour 
of the functional relations for k^u) and k s (u). The known scaling relations [21] aj, = 2 — 2v, 
a s = ab + v are consistent with this relation and can be used to infer the value of the correlation 
length exponent v. These relations have been confirmed explicity for the eight-vertex [jjj and 
the CSOS |H models for which the exponent v is known. For the present models, in regime III 
we thus expect 
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In regime II 
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These results remain to be confirmed via a direct calculation of the correlation length. 

Our results are consistent with a number of partial checks: 
• Regime III 

(i) There is an equivalence at level I between the B^ model and the degree 2 fusion model. 
Formally setting n 



1 in the Bn model we see that the bulk free energy of the level I B^ 



model agrees with the result obtained for the level I degree 2 fusion A\ ' model The bulk 
free energy of the Bn^ model at the critical point (p — > 0) is consistent with the result obtained 



from the string hypothesis [23 
(ii) Formally setting n 

twice those of the A± RSOS model. This is due to the fact that Dip 



2 in the Dn^ model, the Dip bulk and surface free energies agree with 



A [ p @A\^>. The bulk 

free energy of the Dn' model at the critical point is also consistent with the result from the 
string hypothesis. Further we can check that the exponent 2/v of the D^p model is consistent 
with the result from the thermal scaling relation if we identify the dimension of the generalized 
(1, 3) operator in the conformal field theory (2f| as playing the role of the thermal operator. 
• Regime II 
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(i) The bulk free energy of the level / B^ model agrees with the degree 2 fusion model 
with level I. 

(ii) The results of the level I cj 1 ^ model are consistent with those of the (7 + l)-state model. 
The results of the level 1 Cn^ model are also consistent with those of the (n + l)-state A± 
model in regime III with nome p 2 . This latter correlation length exponent has been recently 



obtained directly for the ABF model in regime III [26|. 

(iii) The bulk and surface free energies of the D% model agree with twice those of the 
RSOS model. 
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Table 1 
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